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Abstract Galaxy clusters have their unique advantages for cosmology. Here we collect 
a new sample of 10 lensing galaxy clusters with X-ray observations to constrain cos- 
mological parameters. The redshifts of lensing clusters lie between 0.1 and 0.6, and the 
redshift range of their arcs is from 0.4 to 4.9. These clusters are selected carefully from 
strong gravitational lensing systems which have both X-ray satellite observations and 
optical giant luminous arcs with known redshift. Giant arcs usually appear in the cen- 
tral region of clusters, where mass can be traced with luminosity quite well. Based on 
gravitational lensing theory and cluster mass distribution model we can derive an Hubble 
constant independent ratio between two angular diameter distances. One is the distance 
of lensing source and the other is that between the deflector and the source. Since angu- 
lar diameter distance relies heavily on cosmological geometry, we can use these ratios 
to constrain cosmological models. Meanwhile X-ray gas fractions of galaxy clusters can 
also be a cosmological probe. Because there are a dozen parameters to be fitted, we in- 
troduce a new analytic algorithm, Powell's UOBYQA (Unconstrained Optimization By 
Quadratic Approximation), to accelerate our calculation. Our result proves that this algo- 
rithm is an effective fitting method for such continuous multi-parameter constraint. We 
find an interesting fact that these two approaches are sensitive to J7a and £1 m separately. 
Combining them we can get quite good fitting values of basic cosmological parameters: 
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1 INTRODUCTION 

Cosmic acceleration expans ion was discovered initially from supernovae observations dRiess et al.l 
ll998llPerlmutter et allll999l) . Then the cosmic microwave background (CMB) anisotropy power spec- 



trum (Sp ergel et al.l 120031) confirmed that our universe is nearly fiat and matter densi ty is relative l y low . 



Nowadays s uch conclusions have been supported by more precis e supernova data (IRiess et al. ,_ 2004 



Davis et all 120071: iKowalski et al.ll2008l) , and CMB observations (ISpergel et al.l 120071: iKomatsu et al 



2009b . Furtherm ore, many other ind ependent work, like the light elements abundance from Big Bang 



Nucl eosynthesis ( Buries etail l200ll) . the bar yon acoustic oscill ations(BAO) detected in S PSS sky sur- 



vey ( Eisenstein et al.l 20051) 7 radio galaxies (iDalv et al.L 120091) . and gamma-ray bursts dAmati et all 



2008) etc. also give consistent results. If there is a component with negative pressure, dubbed dark 
energy, full of our universe, such acceleration can be explained within existing theoretical framework. 
Then various models were proposed as candidates of dark energy, such as the typical dynamical scalar 
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field quintessence dCaldwelletallll998l) . phantom correction s dCaldwelil2002h. a joint quint om sce- 
nario dFeng et al. 1 120051) and chaplygin gas (IGorini et al 1 120051: IZhul 12004 IZhang & ZhuLl2006l) etc. On 
the other hand there are still many other astronomers doubting the existence of such a strange mate- 
rial. They are trying to find the other way to understand this accelerat ing universe. Modified Newtonian 
dynam ics(MOND) dMilgromlfeOQll) . Modified Fried mann Equat i on dFreese & Lewisl 120021: IZhu et all 
l2004bl) . Dvali-Gabadadze-Porrati(DGP) mechanism dDvali et al.L l200ol) and so on are all beneficial at- 
tempts. But until now none of these models has overwhelming advantages. We need more observational 
evidence. Besides updating the precision of current data, we also keep trying new probes, for example, 
galaxy clusters. 

Galaxy clusters, as the largest dynamic systems we have known in the universe, are keeping deep im- 
print of the big bang. Its correlation function provides direct measurement of matter d istribution power 
spect rum. Since the BAO peak has already been found with luminous red galaxies dEisenstein et all 
120051) . When there are enough galaxy clusters, we can do the same measurement in a much larger 
scale dBorgani & Guzzol 1200 1[). And their mass distribu tions at different redshifts can be described by 
the Press-Schechter function ( Press & Schechteilll974l) . This relation reflects the linear growth rate of 
density pertur bations. With such co nnection clusters can provide constraints on matter and dark en- 
ergy densities dBorgani et allll999b . What's more, hot gases of galaxy clusters also interact with cos- 
mic microwave background photons and distort its spectrum . Such phenomenon is called the Sunyaev- 
Zel'dovich effect (SZ effect in short) dSunvaev & Zeldovichll 19721) . Combining this effect with the ob- 
servations of corresponding clusters' X-ray luminosity, we can m easure the Hubble constant for a certain 
cosmology or give a rough estimate of cosmological parameters dMason et all 1200 illReese et all 120021 : 
ISchmidt et all |2004 iJones et all 120051: iBonamente et all 120061: lzhu & Fuiimotoll2004 . For more de- 
tailed reviews we suggest Rosati et al.l d2002l) . IVoit ( 2005 ) and Borgani ( 20061) . 

The two methods adopted here are based on the physical structure and properties of individual 
clusters. They can also give a good estimate of cosmological parameters as you will see. The first one 
comes from strong lensing arcs. From X-ray luminosity and temperature we can model mass distribu- 
tion of a cluster's mass halo. Giant arcs generated by high concentrated gravity of a galaxy cluster are 
perfect indicators of its surface mass density. Then we can derive an observational value to constrain 
cosmo logical models. This method was first used bv lSerenol d2002l) . and improved in lSereno & Longol 
d2004l) . We collect a new data set from literature and an online database BAX. Several effective criteria 
are introduced to rule out improper clusters. In the end 10 clusters are selected to ma ke a new sample. 
This n ew set can give an interesting result compared with the last one supplied by ISereno & Longol 
d2004l) . We will discuss it in detail in section 2; The sec ond way is based on the assum ption of con- 
stant x-ray gas mass fraction. This method is developed by I Allen et "aIl d200lLl2004ll2008l) and has been 
proved to be effective. Due to the complex physical mechanism involved, there are many parameters 
to fit. To search multi-parameters space more effectively, we introduce a new optimiz ation algorithm 
UOBYQA (Unconstrained Optimization By Quadratic Approximation) dPowelll 12002) to marginalize 
external parameters. This algorithm has been widely accepted in mathematical field. This should be its 
first application in cosmology. The algorithm is summarized in section 3; The combined analysis and 
discussions are presented in the last section. 



2 LENSING CLUSTER 
2.1 Theoretical foundation 

Gravitational lensin g is one of successf u l predictions of general rela tivity. After it is confirmed by QSO 



observation in 1979 ( Walsh et allll979l) . |Paczvnski & Gorskildl98ll) tried to use these lensing images as 



indicators to estimate cluster mass and constrain cosmological constant. But it is really hard to find the 
deflectors in lensing events o f point-like sources. So this approach was seldom used until recent years 
dFutamase & Yoshidal l200ll). Later on , giant arcs around galaxy clusters were discovered in clusters 
A370 and C12224 dLvnds & Petrosianl Il986l). Thev can also be use d to constrain clusters' projected 



mass, and clusters are easier to observe. Breimer & Sanders! dl992l) did their pioneering work. They 
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estimated the virial mass of cluster A370 with nearly 30 member galaxies' velocities, compared with 
lensing condition and got an original estimate of cosmological parameters. But member galaxies are 
discrete and velocity dispersion of distant objects is usually hard to obtain. Then lSereno & Longold2004l) 
use continuous X-ray luminosity instead. When a galaxy cluster is relaxed enough, the pressure of its 
hot gas c an balance its self-gravity. In t his ca se, we can use hydrostatic isothermal spherical symmetric 
/3-model (ICavaliere & Fusco-Femiano[|l976l) to describe the intracluster medium(ICM) density profile: 



n e (r)=n ea (l + ryrty Wx/2 , (1) 

where n e o is the central electron density, fix describes the slope, and r c stands for the core radius. 
Assuming all the gases have the same emissivity of hot bremsstrahlung radiation, in other words, isother- 
mal (with the temperature T\), the gravity of relaxing cluster and its gas pressure should balance each 
other according to the hydrostatic equilibrium condition. With the approximation of spherical symmetry 
we can estimate mass distribution with gas density, which comes from luminosity fitting result. The 
cluster mass profile can be given in this form: 



MM = 3My/3x r " (2) 
Gfim p r 2 + r 2 ' 

ks is the B oltzmann constant, m p is the proton mass, and /i is the mean molecular weight, which is 
0.6 usually dRosati et all 120021) . Then we can get the proje cted surface mass den sity. Combining the 
result with the critical surface mass density of lensing arcs dSchneider etal.ll 1992b . a Hubbel constant 
independent ratio can be expressed by several observational parameters: 



Dds 



\im v c 1 
obs 6ir ksTxPx 



where Tx, Px, and 8 C are all come from the X-ray data fitting. 8 is a dimensionless angular variable, ra- 
dius divided by angular diameter distance to the cluster (r/Dd). Under flat Friedman- Walker metric,the 
angular diameter distance between an observer at Zd and a source at z s is not equal to D s — Dd- It should 
be integrated from Zd to z s as below : 

Dds = > / -=7^. (4) 

H (l + z s ) J Zd E{z) 

The position of tangential critical curve t is usually deemed equal to observed arc position 9 arc . 
Considering the deviation of extended lensing source position, deflecting angl e has a slight diff erence 
with arc radius angle, 8 t = e9 arc . The correction factor is e = (1/ \A-2) ± 0.04 dOno et all 1 1999b . Then 
the x 2 test can be carried out between observational data and theoretical models. 



2.2 Data selection 



The number of the clusters with arcs is still very limited, and just a small part of these arcs have known 
redshift. At first, we refer to the sample of ISand et al.l (120051) to look for arcs with redshift. Their cat- 
alogue contains 104 tangential arcs from 128 clusters, and only 58 arcs from 27 clusters have redshift 
values. When there are several arcs from the same source (with the same redshift), we prefer to choose 
the farthest one in general. Because strong lensing arcs usually happen in the very central part of galaxy 
clusters, giant arcs are always close to the core of clusters. While the beta model brightness profile usu- 
ally has better fitting results at outskirts. So we choose the farther part to avoid possible substructures 
and decrease the fitting deviation of f3 model. Then arc H5b (26.3 arcsec) of A2390 is adopted instead 
of H5a (20.7 arcsec). For the arcs from different sources around one cluster, we treat them dependently. 
They will all be adopted as long as they can satisfy our criteria given below. 



4 



H. Yu & Z.-H. Zhu 



Second stage, redshifts and temperatures of these galaxy clusters can be searched out directly 
from online databases, such as CDS (The Strasbourg astronomical Data Center) or NED (NASA/IPAC 
Extragalactic Database) with their full name. Here we choose a new database established especially for 
X-ray galaxy clusters - BAXQ- It provides detailed information of clusters. Following referred literature 
given by BAX we can get fitting parameters j3 and 6 C . To minimize systematic error brought by differ- 
ent work and systems, it is better to limit data s ources within a few articles. In this paper, we use the 
fitting result of Cha ndra ( Bonamente et al. , 120061) . For the clusters not presented by that paper, we refer 
to the catalogusfl of lOta & Mitsudal d2004), which is based o n the observations of ROSAT and ASCA 
satellites. The clusters inherited from Sereno & Longold2004 are all updated in this way. 



Table 1 Sample of Lensing Galaxy Clusters with X-ray Observations 



Cluster 


Arc 


Zd 




Oarci. ) 


kTx(keV) 


Px 


OcO 


n(poO)) 


ref 


3C220.1 


Al 


0.61 


1.49 


20.3 


5.56±1.38 


0.84±0.45 


8.1 ±4.2 


8537 


1 


Abell 2390 


H5b 


0.228 


4.05 


26.3 


9.35±0.15 


0.46±0.01 


12.1±1.4 


59463 


1 


Abell 2667 


Al 


0.226 


1.034 


14.7 


6.15±0.61 


0.52±0.01 


13.4±0.8 


65462 


1 


MS 0451.6-0305 


Al 


0.550 


2.91 


31.8 


10.4±0.7 


0.767±0.018 


33.5±1.2 


4177 


2 


MS 1512.4 


cB58 


0.372 


2.72 


5.1 


3.39±0.4 


0.54±0.06 


8.3±2.4 


52805 


1 


MS 2137.3-2353 


A01 


0.313 


1.501 


15.5 


4.96±0.11 


0.6±0.04 


8.3±1.5 


38603 


1 


PKS 0745-191 


A 


0.103 


0.433 


19.2 


7.97±0.28 


0.52±0.01 


16.4±0.9 


271090 


1 


Abell 68 


COc 


0.255 


1.6 


8.0 


10.0±1.1 


0.721±0.035 


49.6±3.6 


17867 


2 


CL0024.0 


E 


0.391 


1.675 


4.0 


4.38±0.27 


0.41±0.03 


11.1±4.1 


31586 


1 


MS 2053.7 


AB 


0.583 


3.146 


15.1 


4.7±0.5 


0.639±0.033 


15.3±1.6 


6324 


2 



Reference: lOta & Mitsuda (2004); 2 Bonamente et al. (2006) 



Then we collect about 20 clusters with all necessary parameters. But not all of them are approxi- 
mately isothermal, spherical symmetric, or even hydrostatic equili brium. We must ch eck them carefully 
with criteria to eliminate the arcs generated by unrelaxed clusters dSmith et al.ll2003l) . 

2.3 Criteria 

A static cluster should have regular morphology both in optical and X-ray band. In optical band they 
usually present spatial symmetry and regular shape. In X-ray band sharp central surface brightness 
and regular elliptical isophotes are both common characters of dynamically relaxed clusters. But such 
description rely on the resolution of equipment greatly. We need more concrete standards. 

Weak gra v itation al lensing has been used to test X-ray mass distribution many years ago. 
iMahdavi et ail d2008l) confirmed th at they were consist ent with each other at inner part. Although he 
used NFW dark matter halo model dNavarro et al.L[l996h . the results make sure that hydrostatic equilib- 
rium is applicable within radius r25oo (which means the mean mass density is 2500 times of the universe 
critical density at that redshift). We can find from Table [T] that all our arcs appear deep inside this re- 
gion. To avoid possible big disagreement of inner mass estimation, we also compare lensing mass with 
theoretical predictions to check our mass model. According to the strong lensing equation of isothermal 
sphere, the projecting mass is: 



AttG D d D ds 




(5) 



1 BAX: http://bax.ast.obs-mip.fr/ , same with the other database like Simbda, NED, etc., only accept clusters' full name like 
RXJ1347.5-1 145. The abbreviations listed here is only for convenience, can not be used to search directly. 

2 Ota's catalogue wasn't contained in the article, but was put on her own website. 
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arc lensing x10 15 

Fig. 1 Relationship between lensing surface mass density(A'/© /Mpc 2 ) and projection sur- 
face density from X-ray model. Blue asterisks come from standard model (f?M = 0.3,Qa = 
0.7);and red circles are calculated under = 0,£Ia = 1. 

And from equation (0, we can derive a theoretical typical surface density, which is: 

Etfc = 3 k BTxPx 1 
2G\im v 6 C Dd 

Mass values from above two equations should be consistent with each other for a spherical relaxed 
cluster as shown in FigQ] Because this test is not model independent, we don't take it as a selecting tool, 
just use it to examine the data with the best fit model. 

In equation distance between the len and the source (integrated from Zd to z s ) should be always 
smaller than that between the arc source and the observer (integrated from to z s ), although the angular 
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D ds /D s 




0.8 







Fig. 2 3-D Hubble diagram distribution of sample galaxy clusters. X axis is the redshift of 
deflectors, Y the arc, and Z is the second criteria Dd s /D s , the reference surface comes from 
our best-fit model below. The intersecting line between fitting surface and the bottom of the 
box is Zd = z s . And the jagged edge is merely caused by programme. 



distance is not monotonically ascending with redshift. So these clusters should satisfy: Dd s /D s 1 j, s < 1 . 
In Figf2]we can see this more clearly. The points above the box indicate corresponding clusters are not 
self-consistent. They need more precise descriptions for cosmological test. Unfortunately this condition 
rules out half of selected lensing arcs. 

And we also discard the arcs whose positions are too far from characteristic radius (8 arc > 3# c ). 
Such cluster has a relatively small core radius and a much bigger arc radius. In this case, even X-ray 
observations can not trace the matter at the region where lensing arcs existing. The extrapolation results 
of the model for inside region is not reliable any more. They may influence constraint result greatly. 
For example, the most luminosity cluster in the sky RXJ1347 is still experiencing merging process 
dAllen et all |2004 lOta et all 120081) . It will cause reall y serious deviation when counted in. In fact it is 
in charge of the high fl\ value (abo ut 1) in the paper of | Sereno & Longold2004l). T his may also explain 
the systematic bias in the articles of lBreimer & Sanders! (1 19921) and ISerenol (120021 A similar situation 
also happens to C arc (38.1 arcsec) of A2390. So only an inner arc H5b is adopted. 

At last we get a sample of 10 clusters listed in TableQ] The redshifts of these clusters rang from 0.1 
to 0.6, and the farthest arc is H5b of A2390 z=4.05. All these arcs appear in very central region. The 
average surface density within the arcs are listed at n column with the unit of the critical universe density 
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Table 2 Values of D ds /D s from Eq.(3) 



Cluster 


D ds /D s 


a(D ds /D s ) 


3C220.1 


0.611 


0.530 


A2390 


0.737 


0.053 


A2667 


0.837 


0.124 


MS0451 


0.785 


0.087 


MS1512 


0.734 


0.330 


MS2137 


0.778 


0.105 


PKS0745 


0.818 


0.065 


A68 


0.982 


0.225 


CL0024 


0.919 


0.430 


MS2053 


0.968 


0.209 



at corresponding redshift. And data for cosmological fitting are shown in Table [2] Errors are calculated 
with error propagation equation. 



3 X-RAY GAS MASS FRACTION 

There is another way to connect physical characters of galaxy clusters to cosmological parameters. As 
mentioned above, in hydrostatic equilibrium the X-ray gas mass distri bution and the total ma ss profile 
which balance it can both be derived from the surface brightness profile dWhite & Frenkl[T99ll) . Because 
of the large scale o f clusters, their matter content can be taken as a fair sample of the whole Universe 
dWhite etal.lll993l) Since the baryonic fraction of clusters are se nsitive to angular d iameter distance. 
ISasakil dl996l) : iPenl (1 1997b first applied this to cosmological test. lAllen et all d200ll) assume that the 
X-ray gas mass fraction within f25oo is a constant. Th e n the a ngular dis t ance of clusters can be der ived 
based on reference cosmological model dAllen et all 12004 IZhu et all l2004al: lAllen et all 120081) . In 
reference cosmology h = Ho/lOOkms^ 1 Mpc -1 ,£Ia = 0.7 and fi m = 0.3, we have 



MCDM 
J qas 



KAjb(z) n b 



1 



d A A CDM {z) 



jmod 

a A 



1.5 



(7) 



where K is a calibration constant, A equals to (#2^o' M /^250o)'^ 7 stands for non-thermal pressure in 
the clusters; b(z) is the depletion factor with the expression of &o(l + a^z); s(z) models the baryonic 
mass fraction in stars, also expressed as sn(l + a s z). Two weak priors are al so needed here: H ubble 
parameter h = 0.72 ± 0.24 and mean baryon density tt b h 2 = 0.0214 ± 0.006 dAllen et all 12008b. 



Obviously, so many parameters aren't easy to calculate with normal meth ods. 



Allen i 



2008h 



l e t al l ( 

used the popular MCMC (Markov Chain Monte Carlo) programme CosmoMC dLewis & BridleL 120021) . 
For such a nonlinear or non-derivable function, stochastic process can generate right distribution, but it 
need a large point set to draw smooth contour lines. And we usually don't need to obtain all parameters' 
bias at the same time. In this case we only care about cosmological parameters. So we can just constrain 
2 parameters each time and marginalize the rest. In actual calculation we assemble the parameters K, 
7, bo, S*o + l into one factor. Marginalizing them together will not affect our final results. Then we 
use grids to generate reference points on parameters' phase space. It is easy to get contour lines with 
them. We just need to search for the most optimistic value in other parameters' space to accomplish 
marginalization. Considering all the physical processes are still continuous, their functions have definite 
slope, we adopt a new analytical algorithm-direct search algorithm. Its main motivation is using nu- 
merical finite difference to approach directional derivatives. Here we use a ready-made program module 
CONDOR (COnstrained, N on-linear, Direct, paralle l Optimization using trust Region method for high- 
computing loa d function^ dFrank & Bersin 1. 12005b . It is based on the Powell's UOBYQA algorithm 
dPowellLl2002t) . 



3 It is available at http://www.applied-mathematics.net 
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With a sample of 42 galaxy clusters jAllen et all l2008h . we obtain the marginalized 68 percent 
confidence limits of £Im = 0.26 ± 0.04 and fl\ = 0.9ig i8> mini mum y 2 = 41.8. Co ntours are shown 
with dash lines in Fig(3 This result is consistent with the results of lAllen et al. I (l2008l) very well, which 
is x 2 = 41.5, Ojvf = 0.27 ± 0.06 and fl\ = 0.86 ± 0.19. Contour lines are showing with dashed lines 
in Fig [3] 

4 RESULT AND DISCUSSION 

Our sample is not big enough to make a precise constraint to the equation of state. So we use a simple 
cosmological model E(z) = -\/£1m(1 + z) 3 + + (1 — — ^a)(1 + z) 2 . Because these two 
methods use different cluster sample, we fit them separately and sum their \ 2 up to get the final fitting 
results. As can be seen in FigfJ) dash lines give the result of X-ray gas mass fraction, and dash dot 
contours come from lensing clusters. The shaded region shows the combining constrain. The small 
circle inside gives the best fitting values: Q, M = 0.261;q o4 and A = 0.82^°;^ at 68% C.L. . These 
results are in agreement with the basic facts we know from other observations. 

When combining results of these two methods, we find that these two constraining methods coming 
from clusters are sensitive to different cosmological parameters separately. Their contour regions are 
nearly orthogonal. The angular diameter distance ratio estimated from strong lensing arcs relies on Ha 
sensitively. The best fitting value is S^a = 0.74^Q'gg. It's better than the result from constraint of gas 
fraction. But it can not give an significant constraint to m- On the other hand, the gas fraction of 
clusters is more effective in constraining il^i- So when we combine these two methods together, unified 
constraint is better than both of them. Such results are still preliminary. Due to the complexity of clusters 
themselves, there is still a lot of work to do to improve these methods. 

Common radial temperature distribution of clusters make it difficult to use isothermal approxima- 
tion. But in the small scale within arcs temperatures won't change steeply even for cool-core cluster. 
The clusters which have substructures or even are experiencing merging can not be discribed sim- 
ply. At least we can rule them out by our criteria. The X-ray observations of our n ew lensing cluster 
sample come from three different satellites: ROSAT, ASCA (|Ota & Mitsudal |2004 and CHANDRA 
dBonamente et al.U2006l) . The differences of equipments may also cause systematical errors which are 
hard to estimate. Considering previous two telescopes are comparatively out of time, fitting results from 
CHANDRA or XMM may give better results. 

These two methods use two different cluster sample sets with some common clusters (A2390, 
MS2137 etc.). The deviation of different observations may also be counted in indirectly. So it's neces- 
sary to select a unified set to share the same clusters and observational parameters. On the other side, the 
mass profile model used by these methods are different. The strong lensing approach uses the /3-model 
to fit luminosity profile of a cluster, while the gas fraction method uses the NFW model to describe dark 
matter halo. As we have seen in the TableQ] lensing arcs usually appear in the central region of clusters, 
where the NFW model can usually give a better fitting result than that of the isothermal sphere model 
dComerford et all l2006t ISchmidt & AllenLl2007h . If we can unify the two methods with the same mass 
profile model, results will be more convincing. 

Compared with the other cosmological ob servations, our clusters sample is really small, and its 
range of redshift is also limited. According to lYamamoto & Futamasd d200ll) . a data set containing 
more than 20 clusters can constrain dark energy equation of state more precisel y. It won't take a long 
time to achieve that g oal. There are many giant arc survey projects proceeding (iGladders et all 12003b 
iHennawi et alll2008l) . The number of new discovered arcs is increasing rapidly. Their redshift measure- 
ment is only a matter of tim e. And next generation X-ray telescope, eg. International X-ray Observatory 
(IXO) (I White et all 1201 Oh. extended ROentgen Survey with an Ima ging Telescope Array (eRosita ) 
dPredehlet all 1201 Ol) and the Wide Field X-ray Telescope (WFXT) (jlylurrav & WFXT Teaml |2010|) . 



will carry out new surveys more precisely in a much larger field. Future observations will definitely 
enlarge our set and make these methods more powerful. 
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Fig. 3 Dash line gives the result of X-ray gas mass fraction. While dash dot line comes from 
lensing cluster. Two "x" simbols give their best fitting values seperately. The shaded regions 
show the 1 ,2, and 3 a confidence regions of combining constrain, corresponding to A^ 2 values 
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